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The ep — + e'n^n reaction was studied in the first and second nucleon resonance regions in the 
0.25 GeV^ < < 0.65 GeV^ range using the CLAS detector at Thomas Jefferson Nationai Ac- 
celerator Facility. For the first time the absolute cross sections were measured covering nearly the 
full angular range in the hadronic center-of-mass frame. The structure functions cttl, cttt and 
the linear combination ctt + eo"L were extracted by fitting the ^-dependence of the measured cross 
sections, and were compared to the MAID and Sato-Lee models. 

PACS numbers: 



INTRODUCTION 



polarization parameter e: 



The structure of the nucleon and its excited states has 
been one of the most extensively studied subjects in nu- 
clear and particle physics for many years. It allows us 
to understand important aspects of the underlying the- 
ory of the strong interaction, QCD, in the confinement 
regime where solutions are very difficult to obtain. Elas- 
tic electron scattering experiments provide information 
on the ground state of the nucleon, while studying the 
evolution of the transition amplitudes from the nu- 
cleon ground state into the excited states provides insight 
into the internal structure of the excited nucleon. Single- 
pion electroproduction is one of the most suitable pro- 
cesses for studying the transitions to states with masses 
below 1.7 GeV because of the large ttN coupling for these 
states m . The detection of two out of three outgoing par- 
ticles is sufficient to achieve a complete measurement of 
the differential cross sections in order to attempt the ex- 
traction of the amplitudes for the individual resonances. 
The kinematic quantities of the ep — > e'n^n reaction is 
shown in Fig. ^ The virtual photon is described by the 
four-momentum transfer Q^, energy transfer v and the 
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y = E,-Ef, (3) 

where Ei and Ef are the initial and final energies of the 
electron and 9e is the electron scattering angle. The mass 
of the hadronic system is given by: 
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where M is the proton mass. The two hadron production 
angles 9 and cj) are defined in the center-of-mass (cm.) 
reference frame, with 9 being the angle between the out- 
going pion and the direction of the three-momentum 
transfer, and (j) being the angle between the electron scat- 
tering plane and the hadron production plane. The un- 
polarized cross section for single-pion electroproduction 
can be written as Q: 
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FIG. 1: (Color online) Kinematic diagram of single-pion elec- 
troproduction. 



where T is the virtual photon flux, and is the vir- 
tual photoproduction cross section. The ctt, cl, cttt 
and aTL structure functions are bilinear combinations of 
the helicity amplitudes, depending only on the variables 
Q^, W and 9. The analysis of the angular distributions 
provides information for extracting the electroproduction 
amplitudes for different resonances. 

The main tree-level Feynman diagrams contributing 
to the ep c'tt+ti process are shown in Fig. |2] The 
s-channel resonance excitation process is represented by 
the diagram in Fig. The hadronic vertex of this pro- 
cess is known from ttN elastic scattering experiments jH . 
Therefore studies of pion electroproduction can yield the 
evolution of the photocoupling amplitudes describing 
the 'y*NN* vertex. For the purpose of studying the exci- 
tation of nucleon resonances, the other diagrams are con- 
sidered as physical background. The largest non-resonant 
contribution to the cross section comes from the f-channel 
pion exchange diagram, shown in Fig. Although this 
process mainly contributes in the forward region due to 
the pion propagator pole, it still accounts for a signifi- 
cant part of the cross section even at large angles. The 
diagrams in Fig. ISb and Fig. [21i correspond to the s- 
channel nucleon pole and i-channel p-meson exchange 
amplitudes. Sophisticated analysis procedures are nec- 
essary to separate the resonant contributions from the 
non-resonant background, and to extract the resonant 
amplitudes for different overlapping excited states. The 
extraction of resonance multipoles is beyond the scope 
of this paper. In this contribution we describe the ex- 
periment and data analysis, and the extraction of fully 
exclusive and differential cross sections, and determina- 
tion of response functions. 

Electroexcitation of a nucleon resonance can be de- 
scribed in terms of three photocoupling amplitudes Ai/2j 




and 5*1/2. The first two are due to the coupling of 
transverse photons with the proton resulting in a com- 
bined helicity /i = ^ or /i = | respectively. The 6*1/2 
amplitude is present due to the possibility of a longitudi- 
nal polarization for virtual photons. Alternatively, pion 
electroproduction can be described using multipole am- 
plitudes Ei±, Mi± and Si±. The Z-index represents the 
orbital angular momentum of the ttN system, and the 
± sign indicates how the nucleon spin is coupled to the 
orbital momentum. For each excited state the helicity 
amplitudes can be expressed in terms of multipole am- 
plitudes and vice versa . 

Quark models predict that the and ratios for 
the P33(1232) are small at low while perturba- 

tive QCD predicts = 1 and-^^ is independent of 
as — > 00 6]. A transition between these two regimes 
is expected at some finite . At low the deviations of 
these ratios from zero can be interpreted as non-spherical 
deformation of the nucleon or the A(1232) Q. Usually 
these ratios for A(1232) are obtained through measure- 
ments in the ir'^p decay channel with an assumption that 
the uncertainty due to the isospin I — \ background is 
negligible. High quality data in the -K^n channel will 
enable us to separate the isospin I — \ and I = \ com- 
ponents of the transition form-factors for the ^33(1282) 
and to determine these ratios with smaller uncertainties 
coming from non-resonant contributions. 

The second resonance region is dominated by the 
three known isospin I — \ states, Pii(1440), 1513(1520) 
and S'ii(1535). These resonances, produced in electron- 
proton scattering, are twice as likely to decay through 
the TT'^n channel than through Tr^p. Therefore, cross sec- 
tion measurements of the ep — > e'vr process are crucial 
for understanding the properties of these states. The na- 
ture of the Pii(1440) resonance is not understood in the 
framework of the constituent quark model (CQM) 0, 
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and there are suggestions that the Roper resonance may 
be a hybrid state j^j or a smaU quark core with a large 
vector meson cloud The evolution of the A^ji 

photocoupling amplitude for the Roper is predicted to be 
different for 3-quark and hybrid states. Previous analyses 
0,01 indicate a rapid fall-off of Axj2 between Q-^ 
and = 0.5 GeV^, therefore high quality data in this 
region will be very valuable in understanding the nature 
of the Pii(1440). 

The experimental data for the A1/2 transition ampli- 
tude for S'ii(1535) show a significantly slower fall-off 
than predicted by constituent quark models. Most of 
these results are obtained through analysis of 77-meson 
electroproduction data, where there can be no / = | 
background. The proximity of the S'ii(1535) mass to the 
ry-production threshold complicates the analysis of the 
data. High quality single tt^ data currently exist only at 
the photoproduction point, and there is very little data 
for non-zero Q^. The results from analyses of pion and 77 
photoproduction data are significantly different 0| , and 
the source of these discrepancies is still not understood. 
New electroproduction data will allow for a similar com- 
parison between the results from the two channels from 
CLAS to check the consistency of the analysis frame- 
works. These data will also allow for a future combined 
analysis of pion and rj production data, which will pro- 
vide more stringent constraints on the fit. 

Until now there have only been three experiments 
[1^ I15L measuring single 7r+ electroproduction cross 
section in the resonance regions in this range of Q^. In 
all of these experiments the lack of angular coverage in 
the center-of-mass reference frame significantly reduced 
the sensitivity to the resonant amplitudes. The aim of 
the present experiment is to provide differential cross sec- 
tions for the TT+n channel over a large kinematic region 
and with high statistical accuracy, that can be used to- 
gether with other channels to obtain more reliable results 
on the resonance photocoupling amplitudes. 



EXPERIMENT 

The measurement was carried out using the CE- 
BAF Large Acceptance Spectrometer (CLAS) 0| at the 
Thomas Jefferson National Accelerator Facility (Jeffer- 
son Lab), located in Newport News, Virginia. CLAS is 
a nearly 47r detector, providing almost complete angular 
coverage for the ep e'n'^n reaction in the center-of- 
mass frame. It is well suited for conducting experiments 
which require detection of two or more particles in the 
final state. Such a detector and the continuous beam pro- 
duced by CEBAF provide excellent conditions for mea- 
suring the ep — > e'TT^n electroproduction cross section by 
detecting the outgoing electron and pion in coincidence. 




FIG. 3: Three dimensional view of CLAS. 



Apparatus 



The main magnetic field of CLAS is provided by six 
superconducting coils, which produce an approximately 
toroidal field in the azimuthal direction around the beam 
axis. The gaps between the cryostats are instrumented 
with six identical detector packages, also referred to here 
as "sectors" , as shown in Fig O Each sector consists 
of three regions (Rl, R2, R3) of Drift Chambers (DC) 
[lij to determine the trajectories of the charged particles, 
Cerenkov Counters (CC) 1 191 for electron identification. 
Scintillator Counters (SC) j2d| for charged particle iden- 
tification using the Time-Of-Flight (TOF) method, and 
Electromagnetic Calorimeters (EC) [2l| used for elec- 
tron identification and detection of neutral particles. The 
liquid-hydrogen target was located in the center of the de- 
tector on the beam axis. To reduce the electromagnetic 
background resulting from M0ller scattering off atomic 
electrons, a second smaller normal-conducting toroidal 
magnet (mini-torus) was placed symmetrically around 
the target. This additional magnetic field prevented the 
M0ller electrons from reaching the detector volume. A 
totally absorbing Faraday cup, located at the very end 
of the beam line, was used to determine the integrated 
beam charge passing through the target. The CLAS de- 
tector can provide ^ < 0.5% momentum resolution, and 
~ 80% of 47r solid-angle coverage. The efficiency of de- 
tection and reconstruction for stable charged particles in 
fiducial regions of CLAS is e > 95%. The combined infor- 
mation from the tracking in the DC and the TOF systems 
allows us to reliably separate protons from positive pions 
for momenta up to 3 GeV. 
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FIG. 4: Energy deposited by the electron candidates in the 
electromagnetic calorimeter versus their momenta. The black 
hnes show the cut applied for the electron identification. 



Data taking and data reduction 

The data were taken in the spring of 1999 as part of 
the experimental program of the CLAS collaboration. 
The CEBAF 1.5 GeV electron beam was incident on 
a 5-cm long liquid hydrogen target at 20.5 K temper- 
ature. The data were taken at 3 nA nominal beam cur- 
rent, with ±0.04 nA current fluctuations, at luminosities 
of ^ 4 X lO"'"^ cm^^s^^. The size of the beam spot at 
the target was ~ 0.2 mm, with position fluctuations of 
±0.04 mm. The main torus current was set at 1500 A, 
which created a magnetic field of about 0.8 Tesla in the 
forward direction. The magnetic field of the spectrom- 
eter is significantly lower at large angles. The CLAS 
event readout was triggered by a coincidence of signals 
from the electromagnetic calorimeter and the Cerenkov 
counters in a single sector, generating an event rate of 
~ 2 kHz. The total number of accumulated triggers at 
these detector settings was about 4.5 x 10^. The raw data 
were written onto a tape silo of the Jefferson Lab Com- 
puter Center. During the off-line processing each file was 
retrieved from the tape silo and analyzed to produce files 
for general use containing 4-vectors of the reconstructed 
particles. ROOT [i^l files, containing the specific in- 
formation relevant for single tt"*" electroproduction, were 
created and stored on a disk. These files were further 
analyzed to extract the differential cross sections for the 
ep — > e'TT^n reaction. 



Particle identification 

One of the key issues in electron scattering experiments 
is the ability of the detector to reliably identify elec- 
trons. Electron identification at the trigger level was ac- 
complished by requiring a minimum amount of deposited 



TABLE I: Cuts for the geometrical matching in the offline 
analysis software. 



energy in the electromagnetic calorimeter in coincidence 
with a signal in the Cerenkov counter in the same sec- 
tor. Additional requirements were applied in the off-line 
analysis to select events containing an electron. First, a 
geometrical matching was required between the EC and 
CC hits and the associated negatively charged tracks in 
the drift chambers. The values of the geometrical cuts in 
the software are given in Tabled A sampling fraction cut 
was imposed on the dependence of the EC visible energy 
on the momentum to reject the background coming from 
negative pions (see Fig. ^ . The electron identification in 
the off-line analysis can be summarized by: 



EID = TRK (g) CC (g) EC (g) SF, 



(8) 



where TRK stands for track reconstruction in the drift 
chambers, CC and EC are the Cerenkov counter and 
the calorimeter geometrically matched hits and SF is 
the sampling fraction cut described above. In order to 
avoid inefficiencies due to the trigger threshold in the 
electromagnetic calorimeter, only events containing an 
electron with at least 500 MeV momentum were used in 
the analysis. In addition, fiducial cuts, discussed later, 
were applied to select only electrons in the regions where 
the Cerenkov counter efficiency was greater than 92%. 
The final cross sections were corrected for the remaining 
inefficiency of the Cerenkov counters 0| . 

Charged hadron identification in the CLAS detector 
is accomplished using the momentum determined from 
the tracking and the timing information from the scin- 
tillation counters. Fig. shows the distribution of posi- 
tively charged particles at 1.515 GeV electron beam en- 
ergy plotted versus velocity f3 and momentum P. Bands 
due to positrons, pions, protons and deuterons can be 
easily identified. At low momentum the muon band is 
visible as well. The deuterons are produced from elec- 
tron scattering on the aluminum windows of the target 
cell. All positive particles in the region outlined by the 
dashed lines were considered as tt"*". Positrons can be 
separated from pions at low momenta, but at higher mo- 
menta the pion and positron bands merge. Background 
due to muons and positrons is significantly reduced by 
the missing mass and vertex cuts described below. The 
remaining contamination is evaluated as a systematic un- 
certainty. 
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Momentum corrections 

When extracting the resonant parameters for excited 
states, it is important to have the correct value for the in- 
variant mass of the hadronic state. Therefore, it is neces- 
sary to measure the electron momentum with high accu- 
racy. For this reason additional corrections were applied 
to the electron momentum reconstructed by the standard 
CLAS software package. These corrections were deter- 
mined using elastic scattering events from the same runs 
that were used in the single pion analysis. It was found 
that the missing mass determined from the elastically 
scattered electrons is typically ~ 5 MeV below the pro- 
ton mass Mp = 0.938 GeV. Assuming that the scattering 
angle of the electron is measured correctly, and using: 



2MEi - {W^ - A'P 
2M + AE^ sin^ % 



the momentum correction factor can be found as: 
^ ~b7 ~ 2ME~' 



Cn 



(9) 



(10) 



where Ei is the electron beam energy and W is the mea- 
sured recoil mass. This quantity was calculated for dif- 
ferent bins in 0^ € [15°, 55°] and </>* e [-30°, -f 30°] in the 
laboratory frame for each sector and stored in a look-up 
table. The azimuthal angle (/)* is defined within a sector, 
with 0* = corresponding to the mid-plane of the sector. 
The momenta of the electrons from single 7r+ production 
data were corrected using this table on an event-by-event 
basis. This procedure relies on the fact that the relative 
momentum offset is independent of W for a fixed value of 
Oe- It was found that after these corrections were applied, 
the neutron peak in the missing mass of the ep —^ e'lr^X 
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FIG. 6: Position of the missing mass peak in the ep e'-K^X 
reaction versus <j)e in the lab frame: a) - before, b) - after 
momentum corrections are applied. The error bars represent 
the width of the distribution around the neutron mass. 



reaction was within ±2.0 MeV of the neutron mass for 
1.1 GeV< W < 1.6 GeV range. Fig.Elshows the differ- 
ence between the missing mass in ep — > e'n'^ X reaction 
and the neutron mass with and without momentum cor- 
rections. The six gaps between the points are due to the 
six coils of the magnet. The dependence of the peak po- 
sition on (pe is practically eliminated by this procedure, 
and the peak is located much closer to the neutron mass. 



Fiducial cuts 

Although CLAS is a nearly Air detector, it still contains 
significant inactive volumes without particle detectors. 
In addition, some of the detection inefficiencies in the ac- 
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FIG. 5: Distribution of the number of positively charged par- 
ticles versus f3 and p. The visible bands are due to positrons, 
muons, pions, protons and deuterons. All positive particles 
within the area outlined by the dashed lines are considered as 
7r"'"'s in this analysis. 
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FIG. 7: Cerenkov counter efficiency versus the 6^ and (/f>e elec- 
tron angles in the laboratory frame for Sector 4. The black 
curves indicate the outer edges of the electron fiducial regions. 
Each momentum bin is 200 MeV wide. 
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tive volumes are not adequately reproduced in the detec- 
tor simulation software. These areas are near the edges of 
the electromagnetic calorimeter, Cerenkov counter mir- 
rors, the main torus and mini-torus coils, regions with 
broken wires in the drift chambers, and malfunctioning 
phototubes in the TOF system. To eliminate events with 
particles traveling through these regions, a set of fidu- 
cial cuts was developed. For electrons the main bound- 
ary of the fiducial region was defined by the efhciency of 
the Cerenkov counters and the edges of the electromag- 
netic calorimeter. Fig. [7| shows the dependence of the 
Cerenkov counter efficiency versus the 9e and (j)e electron 
angles in the laboratory frame for six 200 MeV wide mo- 
mentum bins in Sector 4. Due to the optics design of 
the CLAS Cerenkov counters there are areas with 
relatively lower efficiency shown with the lighter shade. 
These features are difficult to implement in the detec- 
tor simulation. Only events in the regions within the 
black curves and with the Cerenkov counter efficiency 
above 92% were used in the analysis. An additional set 
of geometrical cuts was applied to reject electrons hitting 
malfunctioning scintillator counters or traversing regions 
with missing or inefficient wires in the drift chambers. 

Two sets of fiducial cuts were used to define the outer 
boundary of the fiducial regions for the positive pions. 
The first set, similar to the electron cuts, was defined in 
such a way that the (j)^^^ distributions of the number of 
events be uniform within the fiducial region. The sec- 
ond set of cuts was applied to ensure equivalent solid 
angle coverage for pions in the Monte-Carlo simulation 
and the real data. This mismatch was due to distortions 
of the minitorus coils which were not implemented in the 
detector simulation package. As in the case with the 
electrons, tracks in the regions with malfunctioning scin- 
tillator counters or broken wires were rejected by another 
set of fiducial cuts. 
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FIG. 8: Missing mass spectrum (a) and the distribution of 
the events versus Zt^ — Ze (b). The arrows represent the cuts 
applied in the data analysis. The solid lines are gaussian fits 
to the data. 
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FIG. 9: Distribution of the number of single tt^ events versus: 
and W - left; cj) and 9 center-of-mass angles - right. 



Kinematic cuts 

The exclusive final state was selected by detecting the 
outgoing electron and the 7r+ , and by requiring that the 
missing particle be a neutron. The missing-mass spec- 
trum in Fig. shows a prominent neutron peak as well 
as events from the radiative tail and from multi-pion pro- 
duction channels. The arrows indicate the cuts used in 
the analysis. The number of rejected events in the tails 
is recovered by imposing the same cuts on the simulated 
events in the acceptance calculations. 

GEANT-based Monte-Carlo studies showed that about 
18% of the positive pions decay in-flight into ^'^v^. Most 
of the momentum of the original pion is carried by the 
which is, therefore, often detected and reconstructed 
as a 7r+ with a significantly different momentum vector. 
In order to reduce the number of the events with decaying 
pions, a vertex cut | — Z^, | < 2 cm was applied on the 



difference of the Z-coordinates along the beam-line for 
the electron and the pion tracks in the same event (see 
Fig. IHId) . This led to a reduction in the number of events 
with decaying tt"*" to 4% with less than 1% losses in the 
number of events when the pion did not decay. 

The kinematic coverage of this experiment is shown in 
Fig. 1^ The grating on the figures shows the binning of 
the data. In the first resonance region this experiment 
covers a range from 0.25 GeV^ to 0.65 GeV^, while 
in the second resonance region the upper boundary of 
the coverage is reduced to 0.45 GeV^. The angular 
coverage in the hadronic center-of-mass frame is nearly 
complete, with the exception of the region > 140°. This 
limitation at larger angles is related to the fact that the 
CLAS coverage for charged particles is limited to 140° in 
laboratory frame. The number and the sizes of the cross 
section bins are given in Table ITTl 



8 



Variable 


# of bins 


Lower limit 


Upper limit 


Width 




4 


0.25 GeV" 


0.65 GeV^ 


0.10 GeV^ 


w 


25 


1.1 GeV 


1.6 GeV 


20 MeV 


e 


12 


0° 


180° 


15° 




12 


0° 


360° 


30° 



TABLE IL The number and the sizes of the data bins. Values 
for the limits indicate the upper and lower edges of the bins, 
rather than the bin centers. 



Acceptance corrections 




150 

e (deg) 
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In order to relate the experimental yields to cross sec- 
tions, acceptance correction factors were calculated us- 
ing the Monte-Carlo method. The GEANT-based de- 
tector simulation package GSIM incorporated the survey 
geometry of CLAS, realistic drift chamber and timing 
resolutions along with missing wires and malfunctioning 
photomultiplier tubes. Because CLAS is a complicated 
detector covering almost Att of solid angle, it is virtually 
impossible to separate the efficiency calculations from 
the geometrical acceptance calculations. In this work 
the term acceptance correction refers to a combined cor- 
rection factor due to the geometry of the detector and 
the inefficiencies of the detection and reconstruction. It 
is defined as the ratio of the number of reconstructed 
Monte-Carlo events to the number of simulated events in 
a given bin: 

A^^. (11) 

AT ■ ^ ' 

^ ^ Sim 

With this definition of the acceptance it is desirable to 
have a realistic physics model in the event generator be- 
cause of the finite bin sizes and bin migration effects, 
which are described later. In this work the MAID2000 
model [i^, which reasonably reproduces both p7r° 
and the current n7r+ CLAS data, was used as an input to 
the Monte-Carlo event generator. The simulated 200 mil- 
lion events were processed using the same software pack- 
age and analyzed with the same cuts that were applied to 
the real data. An acceptance table with 8 x 25 x 24 x 48 
bins in Q^, W, and (j>, respectively, was calculated us- 
ing the definition in Eq. (|ll|l . The fine binning of the 
acceptance look-up table reduces the model dependence 
of the cross sections. The statistical errors for the ac- 
ceptance corrections were estimated using the binomial 
distribution: 



' Ai)in ( 1 — ) 



N, 



(12) 



gen 



where Ngen is the number of the Monte-Carlo events gen- 
erated in the bin. These errors are included in the sta- 
tistical error of the final cross sections. 

The acceptance of CLAS for single-pion electroproduc- 
tion at Q2 = 0.3 GeV^ and W = 1.23 GeV is shown 



FIG. 10: Sample plots of acceptance corrections versus 9 
and (j> pion angles in the center-of-mass frame in the = 
0.3 GeV^ and W = 1.23 GeV bin. The ei-dependence (a) is 
shown at = 116.25° and the (^-dependence is at = 108.75°. 
The width of the curves represent the statistical uncertainty 
for the acceptance. 



in Fig. mi The ^-dependence of the acceptance (see 
Fig. llOh .) exhibits a dip at ~ 45°, which is due to the 
forward beam pipe. Six sectors of CLAS can be clearly 
identified in the plot showing acceptance versus (^-angle 
(see Fig. llOb ) . The width of the curves in these graphs 
represent the statistical error bands. Since a single cross 
section bin contains sixteen acceptance bins, the contri- 
bution of the acceptance statistical error to the total un- 
certainty of cross sections is on average approximately 
four times smaller than the errors seen in these plots. 



Radiative corrections 



In addition to processes which result in the exclusive 
e'TT+n final state, there are radiative processes repre- 
sented by Feynman diagrams similar to the original single 
photon exchange diagrams, but with an additional pho- 
ton leg, that also contribute to the cross sections. The ex- 
perimentally measured cross sections must be corrected 
for such processes, also known as internal radiation. The 
radiative cross section for an exclusive process can be 
written as ^3 ■ 

{AirafdQ'^dW'^dni 



dar 



2(47r)752M^2 
v\/Xw 



dQhdv 



(13) 



where S = 2EiMp, dfi* is the differential center-of-mass 
solid angle of the 7r+, w = Mj. — M"^, L^jj^} and W^i, are 
the leptonic and the hadronic tensors respectively, and 

Xw ^ - ml^ - 

+ sin 9t, sin 9^ cos ((/)^ - ) ) . 



(14) 
(15) 
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FIG. 11: Sample plots of the radiative correction factor R at 
= 0.3 GeV^ The dotted lines are the external corrections, 
the dashed lines are internal corrections, and the solid lines 
show the combined radiative correction. 



Here, 6.,^, (p-^, 9k and (p^ are the pion and radiated 
photon's angles in the hadronic center-of-mass reference 
frame. The integral in Eq. H15|l is taken over the photon 
angles and variable v. 

In addition, there is also a nonzero probability that 
in the presence of the electromagnetic field of the atoms 
of the target the electron will emit one or more photons 
before or after interacting with the nucleus of the target 
(external radiation). The probability of emitting a real 
photon of a particular energy is proportional to the path 
length of the electron in the target material. The size 
of the external radiative corrections for these measure- 
ments was significantly smaller than for internal because 
of the small amount of the target material (t = 0.5% of 
radiation length). 

The internal radiative corrections for the cross sections 
were calculated using the ExcluRad program 25] as mul- 
tiplicative correction factors for each data bin. The ex- 
ternal radiative corrections were done using the Mo and 
Tsai formalism |2^. The unradiated structure functions, 
needed as an input for the correction procedure, were 
calculated using a parameterization of the multipole am- 
plitudes using a fit of these CLAS data based on the 
unitary isobar model . The size of the required cor- 
rections varied up to 55%, depending on the kinematics. 
Fig. shows dependences of the radiative corrections 



on the variables W, 9 and at = 0.3 GeV^ 
The dotted line shows the correction due to external ra- 
diation, and the dashed line is the correction factor ob- 
tained using the ExcluRad program. The solid line is the 
combined correction factor calculated as the product of 
the two. Because of the short length of the target, the 
external radiative corrections are much smaller than the 
internal corrections. 



Corrections for binning effects 

Because of the finite detector resolution and finite bin 
size, the measured values of the cross section in the center 
of the data bin can be distorted by up to 10%. The exper- 
imentally measured quantity is the cross section averaged 
over a full data bin, while usually it is more desirable to 
determine the value of the cross section at the center of 
the bin. To account for such distortions, multiplicative 
corrections were introduced as the ratio of the cross sec- 
tion in the center of a bin to the average cross section in 
that bin: 



B ^ 



da 



d(T 

dn* 



(16) 



The averaged cross sections were evaluated using two 
models: the dependence of the cross sections was 
taken from the MAID2000 model "2^, while the cross 
sections at fixed values of were obtained using a uni- 
tary isobar ,211 these data in the first iteration. 

Normalization 

The integrated charge of the electron beam passing 
through the target was measured using the Faraday cup 
located at the end of the Hall B beam line. It generated 
pulses with a frequency proportional to the beam current 
with 10 Hz per 1 nA linear slope parameter. The calibra- 
tion parameters of this device are known with less than 
0.5% uncertainty. The measured charge was corrected 
for the data acquisition live-time, calculated as the ra- 
tio of the counts from two scalers. These scalers were 
connected to a single 100 kHz pulse generator. One of 
them was ungated, while the other one was gated by the 
data acquisition "live" signal. To ensure the quality of 
the analyzed data sample, software cuts were imposed on 
the live-time, elastic scattering and single tt^ electropro- 
duction rates. The portions of the runs for which these 
quantities were outside of the imposed limits were ex- 
cluded from the analysis, with the corresponding beam 
charge being subtracted from the total charge. As was 
mentioned above, the Cerenkov counter efficiencies were 
parameterized during the calibration procedure (loj , and 
the appropriate corrections were applied to the cross sec- 
tions. The comparison of the elastic scattering cross sec- 
tions versus 9e from CLAS and the model calculation 
using a parameterization for the elastic form factors 
is shown in Fig. 1121 The model cross section includes ra- 
diative effects, according to the Mo and Tsai formalism 
[2^ . The error bars on the data points represent statisti- 
cal uncertainties only. The solid line at i? = 1.015 shows 
the result of fitting a constant to the ratio of the mea- 
sured cross sections to the parameterization l^^. The 
fluctuations around this line can be used to estimate the 
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FIG. 12: The ratio of the measured elastic cross section to 
the parameterization of the world data [2^. The error bars 
represent the statistical uncertainty only. The solid line is 
from the fit of the data points to a constant. 



systematic uncertainty of electron detection and recon- 
struction. 

The contributions from scattering off the target cell 
walls was estimated to be 1.5% using empty target runs. 
This correction factor was applied to the cross section in 
every data bin. 



SYSTEMATIC ERRORS STUDIES 

A number of studies were carried out to estimate the 
systematic uncertainties on the measured cross sections. 
The primary method used in these studies was to vary 
different independent parameters of the analysis to de- 
termine the corresponding change in the resulting cross 
sections and the structure functions. 

Because of the finite bin size, the result of averaging the 
acceptance over an acceptance bin depends on the distri- 
bution of events in that bin. If the physics model used 
in the Monte-Carlo simulation differs from the real data, 
then the averaging over a bin may result in an incorrect 
cross section. The introduced error depends on the shape 
of the acceptance function and the cross sections as well 
as on the acceptance bin size. In addition, because of the 
finite detector resolution, some of the events produced 
in one acceptance bin will be reconstructed in a different 
bin. This may cause significant distortions in the final 
cross section distributions. In order to correctly account 
for these effects, a realistic physics generator and detec- 
tor simulation are required. To estimate the errors of 
the final results due to the model used in the acceptance 
calculations, we calculated the acceptance table with two 
different models. The comparison of the results with the 
two acceptance corrections allowed us to estimate the sys- 
tematic errors due to the physics model in the acceptance 



calculations. 

As was mentioned above, we use missing mass and ver- 
tex cuts to select the single-pion production events and to 
reduce the number of events with decaying pions. These 
cuts cause losses of some single-pion events as well. The 
true number of events is expected to be recovered by ap- 
plying the acceptance corrections by using exactly the 
same cuts on the Monte-Carlo data. The remaining sys- 
tematic errors associated with these cuts were estimated 
by varying the sizes of the windows. The absolute value 
of the cross section variations calculated with different 
cut windows, averaged over (j> at fixed Q^, W and 9, was 
considered as the systematic uncertainty for all (f> for that 
fixed and 61. 

One of the possible sources of systematic errors in this 
experiment is the uncertainty in the normalization. This 
can arise from miscalibrations of the Faraday cup, tar- 
get density variations, errors in determining the target 
length and its temperature along with data acquisition 
live-time and other factors. However, the presence of the 
elastic events in the data set allows us to account for 
the normalization uncertainties of the cross sections by 
comparing the elastic cross sections to the parameteriza- 
tion of the world data j2^. This way we were able to 
combine the normalization, electron detection, electron 
tracking and electron identification errors into one global 
uncertainty factor. A comparison of the measured elas- 
tic cross sections for different 6e and (j>e electron angles 
allowed us to assign a 5.2% global uncertainty due to the 
normalization and electron efficiency uncertainties. 

The systematic uncertainty due to the model used in 
the radiative corrections was estimated by performing a 
second iteration. The radiatively corrected experimental 
cross sections from CLAS were fitted once more, and us- 
ing the fit the new correction factors were calculated and 
compared with the previous iteration. The comparison 
indicated an uncertainty on the order of 2% due to the 
model dependence of the radiative corrections. 

Using the kinematically over-determined reaction ep 
e'Tr~^iT~p allows us to determine the tt"*" eflticiency by de- 
tecting the outgoing electron, tt~ and proton. The effi- 
ciency of the 7T~^ detection can be found as the ratio of the 
number of events where the tt"*" was detected to the num- 
ber of events where the tt"*" was expected to be detected. 
A comparison of the pion efficiency calculated from the 
real data with the efficiency from GEANT-based Monte- 
Carlo simulation lead to a systematic error estimate of 
2.5%. 

In order to estimate the background coming from 
two-pion production, a sample of two-pion Monte-Carlo 
events was processed as if it were the actual data sample. 
The analysis of these events showed that this background 
would contribute less than 1% uncertainty to the differ- 
ential cross sections. The systematic error due to the 
7r+ misidentification was estimated to be about ^ 0.5% 
by varying the cut in the proton-pion separation in the 
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analysis. 

The total systematic error in each bin was calculated 
as the square root of the sum of the squares of these 
different contributions. The size of the systematic errors 
is typically slightly larger than the size of the statistical 
uncertainties and is shown in Fig. ll3l as the shaded bands. 



RESULTS 

Cross sections 

The experimental differential cross section for each 
data bin was determined using the following formula: 



dEfdVLedni 



R.B j: 



1 



ALec 



AQ^AWsm9A9A(j) d{Ef, cos 0e)' 



d{Ef,COS 9e) 



2ME,Ef 
W ' 



(17) 
(18) 
(19) 



where the sum in Eq. 1171 runs over the ctt^ (n) events 
reconstructed in the fiducial region of CLAS. Here A is 
the acceptance correction factor for an event, L is the 
integrated luminosity, Na is Avogadro's number, p is the 
target density, Lt is the target length, Qtot is the in- 
tegrated charge corrected for the data acquisition live- 
time, e is the electron charge, Ecc is the Cerenkov effi- 
ciency correction factor, AQ^, AW, AO, Acj) are the bin 
aiE^'^sl ) Jacobian between the {W, Q^) and 

{Ef, cosOe) sets of variables, and R and B are the ra- 
diative and binning correction factors, respectively. The 
values of all kinematic variables are calculated for each 
particular event, as opposed to being taken at the cen- 
ter of the bin. The virtual photoproduction cross section 
can be obtained, according to Eq. (|SJ|, by dividing the 
left-hand side of Eq. (|17|l by the virtual photon flux F 
factor defined in Eq. Sample plots of the differen- 
tial cross sections compared with models are shown in 
Fig. El The solid line shows the cross sections calcu- 
lated using MAID2003 2^ multipoles with I < 5 (here 
simply referred to as MAID2003 model). The MAID 
model uses effective Langrangian approach to calculate 
the Born background, including lu and p meson calcu- 
lations. The background is unitarized in the K-matrix 
approximation. The resonant amplitudes are determined 
by fitting the world pion production data. The dashed 
line in Fig. 1131 corresponds to the model by Sato and Lee 
j3fll |. This model obtains an effective Hamiltonian from 
the interaction Lagrangian using the method of unitary 
transformations. Due to complicated calculations Sato- 
Lee model only includes the P33(1232) resonance and, 
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FIG. 13: 



Sample plots of n 

l2 



virtual photoproduction cross 
sections at = 0.3 GeV^ for different kinematics. The 
shaded bands represent the systematic uncertainties. The 
solid curve is MAID2003 '2^ , and the dashed curve is a 
calculation from Sato- Lee |3n|. 



therefore, it's validity domain is limited to the first res- 
onance region. The shaded areas at the bottom repre- 
sent the estimated systematic uncertainties. Typically 
the systematic uncertainties are slightly larger than the 
statistical errors. The data and the models are globally 
in qualitative agreement, while in certain regions there 
are quantitative discrepancies. Due to the large number 
of data points, it is more convenient in this paper to dis- 
cuss the structure functions rather than the cross sections 
themselves. The values of the measured cross sections 
will be available from the CLAS physics database [U or 
upon request [s^ ■ 



Structure functions 

The structure functions cttt, '^tl and the linear 
combination ax + ec^ were obtained by fitting the (jj- 
dependence of the cross section to a function of the form: 



F{(j)) ^ A + Bcos(t> + Ccos2(t>. 



(20) 



The large angular coverage of CLAS in the center-of- 
mass reference frame allowed us to extract the structure 
functions up to 145° in the centcr-of-mass angle. The 
W and 9 dependencies of the structure functions are 
shown in Fig. 1141 and Appendix A. The solid curves in 
Fig. [Hand Appendix A are from MAID2003 jl^ calcu- 
lations, while the dashed curves are from the Sato-Lee 
model 30] . The table of the structure functions are pre- 
sented in Appendix B. The error column in the table 
shows the statistical and systematic uncertainties added 
in quadrature. 
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A typical PF-distribution of the ut + ecr^ and ott 
terms features a distinct A(1232) peak, followed by 
the less prominent second resonance region. The 9- 
dependence of ctt + cctl is mostly flat at low values of 
W near the single-pion production threshold. This is 
consistent with i?o+ dominance at low energies. With 
increasing invariant mass, structures are developing at 
6 ~ 90°, which is characteristic of resonance produc- 
tion. Above the first resonance region the 6'-dependence 
of (Tt + ecTi becomes monotonically falling, consistent 
with the t-channel pion exchange mechanism dominance. 

From the plots in Figs. 19,21,23,25 in Appendix A, 
one can see that in the first resonance region the mea- 
surements of (Jtt agree very well with the Sato-Lee [lOl 
model. For these values of hadronic center-of-mass en- 
ergy this structure function is dominated by the A(1232) 
resonance contributions. But the model predictions for 
the (Tt ecTL linear combination disagree with the mea- 
sured values, although the discrepancy is within the sys- 
tematic errors. The present experiment does not sep- 
arate the longitudinal and transverse structure func- 
tions. Given the successful description of the gtt term, 
which can be expressed in terms of helicity amplitudes 
as (Ttt ~ Rg{Hj,H2 — H4H1), one may assume that the 
structure function ctt ^ \Hi\'^ + |i/2p + li^aP + |i?4p, 
which depends on the same helicity amplitudes, may also 
be described reasonably well. Then the discrepancies in 
the sum ctt -t- ecL could be due to incomplete knowledge 
of the non-resonant physical background contributing di- 
rectly to (7l. Therefore, these data can be used to im- 
prove our understanding of the non-resonant background 
in the first resonance region. The measured cttl struc- 
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FIG. 14: Sample plots of structure functions from CLAS at 
= 0.4 GeV^ The solid curves are from MAID2003 j2al 
and the dashed curves are from the Sato-Lee calculations [30| . 
The shadowed areas show the systematic uncertainty. 



ture function is the smallest, and the relative systematic 
uncertainties are large. The predictions for cttl from 
Sato-Lee [3(1 | are in agreement with the measured values 
within the error bars. 

The MAID2003 model, which is a fit to predominantly 
TT^p channel, describes our data surprisingly well, with 
the curves in the plots following most of the features of 
the experimental data. But the absolute values for the 
(Ttt and aT+^cTL structure functions are typically overes- 
timated, especially in the second resonance region. This 
may be indicative of our relatively poor knowledge of the 
Z?i3(1520) and S'ii(1535) strength in the ep e'7r+n 
channel. We also observe a distinct structure in the 9 
dependence of the cttl ampUtude for W > 1.32 GeV, 
which is not reproduced by MAID2003, where the sign 
is in fact opposite (see Fig. I14|l . Inclusion of these data 
in the MAID fit can improve our knowledge of the reso- 
nance parameters, background terms and the branching 
ratios for the states in the second resonance region. 

SUMMARY 

In conclusion, for the first time we have measured 
the unpolarized electroproduction cross sections for the 
ep — !■ e'TT+n process covering a large angular range in the 
center-of-mass frame, and we have extracted the cttt, 
(Ttl and ctt -I- sctl linear combinations of the structure 
functions. The combined statistical and systematic er- 
rors are of the order of 10% in most of the measured 
kinematic region. In the first resonance region the mea- 
sured cross sections and the structure functions are in 
qualitativ e ag reement with the MAID2003 jl^ and the 
Sato-Lee 30] models, with a quantitative discrepancy 
with the MAID2003 model jl^ . In the second resonance 
region MAID2003 overestimates the height of the reso- 
nance peak. Together with pn^ channel these data will 
provide the basis for the analysis of resonance transition 
form-factors in a coupled-channel analysis. 
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APPENDIX A: PLOTS OF THE STRUCTURE 
FUNCTIONS 
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FIG. 15: (Color online) Structure functions versus W at = 0.3 GeV^ Solid curves represent MAID2003 calculations, while 
the dashed curves show the predictions of the Sato-Lee model. Shaded areas represent the systematic uncertainties. 
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FIG. 16: (Color online) Structure functions versus WatQ^ = 0.4 GeV^ Solid curves represent MAID2003 calculations, while 
the dashed curves show the predictions of the Sato-Lee model. Shaded areas represent the systematic uncertainties. 
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FIG. 17: (Color online) Structure functions versus W at = 0.5 GeV^ Solid curves represent MAID2003 calculations, while 
the dashed curves show the predictions of the Sato-Lee model. Shaded areas represent the systematic uncertainties. 
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FIG. 18: (Color online) Structure functions versus WatQ^ = 0.6 GeV^ Solid curves represent MAID2003 calculations, while 
the dashed curves show the predictions of the Sato-Lee model. Shaded areas represent the systematic uncertainties. 
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FIG. 19: (Color online) Structure functions versus cm. e at = 0.3 GeV^ Solid curves represent MAID2003 calculations, 
while the dashed curves show the predictions of the Sato-Lee model. Shaded areas represent the systematic uncertainties. 
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FIG. 20: (Color online) Structure functions versus cm. 6 aX = 0.3 GeV^. Solid curves represent MAID2003 calculations, 
while the dashed curves show the predictions of the Sato-Lee model. Shaded areas represent the systematic uncertainties. 
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FIG. 21: (Color online) Structure functions versus cm. e at = 0.4 GeV^ Solid curves represent MAID2003 calculations, 
while the dashed curves show the predictions of the Sato-Lee model. Shaded areas represent the systematic uncertainties. 
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FIG. 22: (Color online) Structure functions versus cm. 6 aX = 0.4 GeV^. Solid curves represent MAID2003 calculations, 
while the dashed curves show the predictions of the Sato-Lee model. Shaded areas represent the systematic uncertainties. 
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FIG. 23: (Color online) Structure functions versus cm. e at Q2 = 0.5 GeV^ Solid curves represent MAID2003 calculations, 
while the dashed curves show the predictions of the Sato-Lee model. Shaded areas represent the systematic uncertainties. 
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FIG. 25: (Color online) Structure functions versus cm. e at = 0.6 GeV^ Solid curves represent MAID2003 calculations, 
while the dashed curves show the predictions of the Sato-Lee model. Shaded areas represent the systematic uncertainties. 
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